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Abstract 

We study the photoluminescence from a near-surface quantum well in the regime 
of ambipolar tunneling to the surface states. Under steady-state excitation an 
electric field develops self-consistently due to the condition of equal tunneling 
currents for electrons and holes. The field induces a Stark shift of the photolumi- 
nescence signal which compares well with experimental data from near-surface 
GaAs/AlGaAs single quantum wells. 
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For a quantum well built in proximity of an unpassivated surface, tunneling to surface states can 
be a nonradiative recombination channel competitive with photoluminescence. The importance of 
this effect in determining the emission efficiency has been demonstrated experimentally in various 
papers PJfl. Recently we have proposed a quantitative model based on ambipolar tunneling of 
electrons and holes which is applicable to many-well systems in the bulk or single wells coupled to 
surface states ||. In steady-state situations the ambipolar regime, with equal tunneling currents for 
electrons and holes, imposes an electric field to develop pHl. The field induces a peak shift of the 



excitonic recombination via the quantum confined Stark effect [pLOU - Here we specialize the discussion 
to the case of a quantum well coupled to surface states and we compare the theoretical results with 
photoluminescence experimental data in GaAs/AlGaAs material. 

We will use the following notation. The width of the quantum well is a and the width of the surface 
barrier is b. The bottom of the el and hhl bands of the well are E e i and E^hi and G is the generation 
current density of electron-hole pairs in the well. We assume that no pairs are generated within the 
barrier or at the surface. The pairs generated into the well relax almost instantaneously, compared 
to the other relevant time scales, to the lowest band of the well. Electron-hole interaction leads to 
exciton formation. Tunneling from the well to the surface states is due to free electrons and holes only 
0. On the other hand, photoluminescence is restricted only to excitonic recombination in the well. 
If n w and p w are the steady-state concentrations (number of particles per unit area) of electrons and 
holes in the well, the following rate equations hold: 

= G — J e — \n w p w (la) 

= G - J h - Xn w p w (lb) 

= \n w p w - I . (lc) 
The bimolecular generation rate of excitons is assumed proportional to the electron and hole con- 



centrations |TT] . The photoluminescence current density / is proportional to the exciton concentration 
in the well. Transfer of electrons (holes) from the well to the surface states is realized in a non-coherent 
two-step process. Quantum coherent tunneling of electrons (holes) from an occupied state of the el 
(hhl) band of the well to an equal-energy empty state at the surface is followed by relaxation toward 
the lower energy states. When the barrier width b is not too small, the phonon relaxation process 
at the surface is much faster than the tunneling process (current densities J e and Jh) and can be 
neglected. 

The tunneling current densities are approximately proportional to the charge concentrations in 
the well and the proportionality factor, namely the tunneling probability, is generally quite different 
for electrons and holes. Therefore, in a steady-state situation when J e = Jh, the concentrations of 
electrons and holes in the well must be different. The resulting electric field, in turn, affects the 
electron and hole tunneling probabilities. 

Since the tunneling rate depends both on the effective mass of the carriers and the density of states 
at the surface, two cases are possible. When the density of states of the donor-like band at the surface 
is not sufficiently smaller than that of the acceptor-like band, the electron tunneling rate is larger 
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than the hole tunneling rate. In this case electrons accumulate at the surface and in a steady-state 
situation p w > n w . The electric field is directed from the well to the surface and its value is given by 

- m (2) 



EoE r 

where n = p w — n w and e r is the permittivity of the barrier material. A reversed situation, however, 
may happen when the effective mass difference between electrons and holes is overcompensated by the 
difference in the surface densities. 

For a given value of the electric field F, in the first order perturbation theory the tunneling current 
densities are J e = n w /r e and Jh = p w / T h where 
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m\Ve\^)\ 2 AuMf{-^) (3) 



^i(«)fA,(o)/(i). (4) 

z/ e (e) and Vh{ e ) are the densities of states in the donor-like and acceptor-like surface bands respectively. 
Energies are measured from the bottom of the el and hhl bands. The electron and hole Fermi 
energies e e F and e F which appear in the Fermi function / are related to the respective electron and 
hole concentrations at the surface, n s and p s , as explained in the following. Finally, A is the relevant 
transverse area and the matrix elements are evaluated in the Appendix. 

When the electric field and the tunneling rates are known, the solution of the rate equations is 

n W(i + 2i^) 2+ i-(i + 2^) (5a) 

p w = n w + n . (5b) 

Moreover, n s = n and p s = if F > and n s = and p s = —n if F < 0. 

This result allows one to find the steady-state values of the electric field, of the charge concentra- 
tions and of the tunneling rates by a recursive method. Starting from some trial value, the electric 
field, i.e. the charge concentration n, is changed until the condition J e = J h is reached. At this point 
the luminescence current density / is obtained from the equilibrium concentrations of electrons and 
holes in the well. 

Carrying out explicitly the calculations implies the knowledge of the energy distribution of the 
surface states. Inversely we can try to get information on the surface states by fitting experimental 
photoluminescence data. We concentrate on the specific example of an Alo.3Gao.7As surface with a 
nearby GaAs quantum well fHH. 

At energy close to the bottom of the el band of the well the Alo.3Gao.7As surface has only donor-like 
states belonging to the exponentially vanishing Urbach tail 

, e (e) = - F exp . (6) 
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Note that energy is measured from the bottom of the el band. Such states are assumed to be nodal 
hydrogenic wavefunctions [T2| with radius r e fixed by their depth into the gap. Their explicit expression 



is given in the Appendix. We assume that at the top of the gap the state density is the two dimensional 
density of free Alo.3Gao.7As electrons with effective mass m s e . The parameter e e will be considered as 
a fitting parameter. According to Eq. (||) the Fermi energy for the donor-like surface band containing 
n s electrons is 

/ irh 2 n q \ , , 

e e F = AE c + eFb-E el + e e \n — ^ . 7 

V m e£e J 

On the other hand, at energy close to the bottom of the hhl band of the well the Alo.3Gao.7As 



surface has a very high concentration of acceptor-like defect states p3| . We schematize them again by 
nodal hydrogenic wavefunctions |L2] but with radius to be considered as a second fitting parameter. 
These states are assumed to be distributed in energy with constant density v h over an interval AE h 
into the gap. The Fermi energy for the acceptor-like surface band is then 

4 = Ps/vh + A£„ - eFb - E hhl - AE h . (8) 

Due to the high ratio between the hole and the electron surface-state density [|13|] holes accumulate 
at the surface and in a steady-state situation we have F < 0. 

Experimental photoluminescence data are available (see Ref. Q for details) for a well width a = 60 
A , a temperature T = 4.2 K, a photon-pump energy his = 1.608 eV and with an incident power density 
Pi = 0.5 Wcm~ 2 . The absorption efficiency is estimated to be 1%, so we take G = 0.01 Pi/hu. The 
relevant material parameters are |T3J]: AE C = 0.3 eV, AE V = 0.128 eV, m s e = 0.091 m, m™ = 0.067 m, 



= 0.34 m, m being the free electron mass, and e r = 12. Moreover we put A = 6 cm 2 s _1 [[LI]]. We 
assume an acceptor-like surface state density u h = 10 14 cm~ 2 eV _1 |T3J with AEh — 0.5 eV (the results 
we found do not depend crucially on this particular value). 

The free parameters, e e and r^, are fixed by fitting the normalized photoluminescence intensity 
I /loo to the experimental data [1|] obtained for different values of the barrier width b (the normalization 
factor Too is the photoluminescence current density for b — > 00). A least-square-error procedure gives 
the unique solution e e = 12 meV and = 11 A. In Fig. 1 we compare the ratio I / loo, calculated with 
these values, with the experimental data. The agreement is excellent. 

In Fig. 2 we show the electric field value calculated in the situation of Fig. 1 as a function of 
the barrier thickness. In the same figure we show also the field obtained with different values of the 
incident power density Pi, all the other parameters being fixed. It is seen that, for high levels of 
excitation, the field approaches values of order 10 5 V cm -1 and keeps increasing when the barrier 
becomes thinner. 

An important check of the validity of our model is given by the analysis of the Stark shift. The 
self-consistently estimated electric field induces a band bending which modifies the single particle 
levels el and hhl and, therefore, the exciton recombination energy AE P . The Stark shift calculated as 
a sum of the shifts of levels el and hhl is shown in Fig. 3 as a function of the QW excitation (details 
of the calculation will be reported elsewhere). In Fig. 3 we show also the corresponding measured 
energy shifts (dots). The agreement is good in the region (Pj ~ 0.5) where the fitting parameters were 
fixed as explained above and is fairly satisfactory over three orders of magnitude. An improvement 
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should be possible if the surface-state spectrum were known a priori. This is a confirmation that the 
ambipolar tunneling approach provides a reasonably accurate description of the loss of efficiency in 
near-surface quantum wells. 

Appendix 

The surface is defined by the plane z = and the well is in b < z < b + a. Firstly, we consider 
the case of electrons. We assume a rectangular potential profile with left and right discontinuities 
Vi = AE C + eFb/ 2 and V r = AE C for the well where F is the electric field in the barrier region 
< z < b. The electron wavef unction at energy e = from the bottom of the el band of the well is 
given by 

p ( sin S exp {k\[z — b]) z < b 

$™ =0 = —= I sm(k(z -b)+6) b<z<b + a (9) 

^ [ sin(A;a + 5) exp (— k [z — b — a]) z > b + a 



where tiki = ^J2mf(AE c + eFb/2 - E el ), hk = ^J2m™E el and hk = ^j2mf(AE c - E el ). The phase 
shift is 5 = tan _1 (A;/A; 1 ) and the energy is determined by solving 

ka = n- sin' 1 ( hk =) - sin' 1 | —^= ) (10) 
\^2mfVi) \^mfV r J 

The constant C is fixed by normalization 

sin 2 5 sm 2 (ka + 5) a 



C 



2k x 2k 2 

-1/2 



sin(2(A;a + 5)) - sin(2£) \ 1 

4k J ( ' 



The donor-like surface state $^ =0 is approximated by a truncated 2p hydrogenic wavefunction fT2j 

z 



where r = y/x 2 + y' 2 + z 2 . The state is at energy h 2 / (&m s e r 2 ) below the bottom of the conduction band 
for the barrier material where the electron effective mass is m s e . By imposing the condition that this 
energy corresponds to e = 0, we determine the radius r e 

h (13) 



8m s e (AE C - E el + eFb) 



Assuming that the perturbation potential V e is of the order of the conduction band offset, the 
tunneling matrix element between the well and surface states at e = can be evaluated analytically 
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(n\v e \^) = ^ r3J2AI h Csin5 le-^ 
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4{2hr e - If 



32(2hr e - l) 3 
1 + b/r e 



+ e 



-b/2r e 



{b/r e f + 2b/r e 
2k 1 r e - 1 



+ 



4(2k ire - l) 2 32(2k ire - 1 



(14) 



In the case of holes we have a completely analogous situation where the relevant band in the well 
is hhl instead of el and the acceptor-like surface state is given by Eq. (pT2|) with r e — ■> r^. Equation 



(PI) gives the tunneling matrix element for holes with the substitutions r e 
eFb -> -eF6. 
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FIGURES 



FIG. 1. Normalized photoluminescence ratio I /loo of a near-surface well vs the surface-barrier thickness b. 
Dots: experimental data from Ref. [1]; solid line: best fitting in terms of the self-consistent model. Incident 
power density is Pj = 0.5 W cm -2 . 

FIG. 2. Calculated electric field F across the surface barrier vs the surface barrier thickness b for different 
incident power densities Pj. 

FIG. 3. Comparison between the Stark shift of the photoluminescence signal calculated from the model 
(solid line) and measured (dots) vs the incident power density P«. The surface-barrier thickness is b = 80A. 
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